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SUMMARY 



A method has "been developed for making photo elastic 
analyses of three-dimensional stress systems "by utilizing 
the polarization phenomena associated with the scattering 
of light. By this method, the maximum shear and the di- 
rections of the three principal stresses at any point 
within a model can "be determined, and the t-.ro principal 
stresses at a free-hounding surface can he separately eval- 
uated. Polarized light is projected into the model through 
a slit so that it illuminates a plane section. The light 
is continuously analyzed along its path "by scattering and 
the state of stress in the illuminated section is ohtained. 
By means of a series of such sections, the entire stress 
field may he explored. The method was used to analyze the 
stress system of a simple beam in "bending. The results 
were found to he in good agreement with those expected from 
elementary theo ry # 

INTRODUCTION 



A new method for the photoelastic analysis of three- 
dimensional (spatial) stress systems was reported in de- 
tail at the Eastern Photoelastic Conference at Cornell 
University in May 1939 • A "brief description of the method 
is outlined in reference 1. Following the announcement of 
the method, the N.A.C.A. constructed a single polariscope 
suitable for preliminary work and began the study of the 
application of the method to aircraft structures, power 
plants, and other related problems. 

Several investigators have attempted to develop a 
truly three-dimensional method of photoelastic analysis. 
With the exception of the fixation method reported in ref- 
erence 2, these attempts have lacked general usefulness. 
In the present paper, the "basic principles underlying the 
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proposed, throe-dimension.nl method are discussed and a con- 
crete example is worked out to show its application. 



THE PKOTOELASTIC METHOD 



The ohotoelastic method (references 3, 4, and 5) is 
"based on the fact that many materials, notably certain 
plastics, "become doubly refractive ( bi r ofringent ) when sub 
jected to stress. The amount of this birefringence may be 
studied with a polariscooe and the character! sti C8 of the 
stress system obtained. A model of a machine part or of a 
structural part may therefore be made from such a material 
and subjected to an appropriate loading system; and the 
state of stress within it can then be determined. If the 
materials from which the model and its prototype are made 
possess similar (proportional) elastic properties, the 
stress systems will coincide in direction and be propor- 
tional in magnitude provided that the deflections are such 
that the geometry remains substantially undisturbed. 

It has been common in the past to study stresses in 
this manner, but such study has been limited in practice 
to plane stress systems. Models were cut from plates and 
loaded in a single plane* Many problems have been solved 
in this manner and valuable data obtained. 

Models have been made from glass, celluloid, gelatin, 
Bakelite, M.arblette, and rubber, but the material known as 
Bakelite 3T-S1-393 has been generally accepted as the most 
suitable for a wide variety of tests. It remains linearly 
elastic over a wide stress range (0 to 8,000 lb. per sq.in 
its elastic modulus is not too low (600,000 lb. per sq.in. 
and its stress-optical coefficient is high (35 lb. per sq. 
in. per fringe for green light). A discussion of the prop 
erties of photoelastic materials is giv^n in reference 6. 

SPATIAL PHOTOELASTIC ANALYSIS 



The theory of photoelastic analysis as applied to 
three-dimensional stress systems will now be considered in 
some detail. The difficulty in t hr e ^-diraen si onal photo- 
elastic analysis by conventional methods has been that, 
when a beam of light passes through a transparent model in 
which a spatial stress system exists, the relative retarda 
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tion observed "between the components of the emergent beam 
is an integrated effect along the total path traversed. 
No means have hitherto "been available for conveniently dis- 
covering the way in which the relative retardation takes 
place along the path. 

When light is scattered within a material medium at 
right angles to the incident "beam, the scattered light 
(reference 7) is plane polarized. This scattering proper- 
ty may "be used either as the polarizer or the analyzer in 
a polariscope. In Bakelite, this phenomenon is perhaps 
due to the micelles discussed "by let lay i (reference 2 ) . 

In the investigation of the variation in the effect 
of the birefringence on a light beam as it passes through 
a model, a source of polarized light within the model capa- 
ble of motion from point to point might be desirable. In. 
order to produce such a source, an unpolarized beam may be 
projected into the model and the scattered li^ht may be 
viewed from a direction 90° to the beam. This arrangement 
may most conveniently be set up by collimating the beam 
with a lens and subsequently passing it through a slit, 
say one-eighth of an inch wide and as long as desired. 
When this beam is passed through the model, it illuminates 
a plane section. Viewed along any path normal to the illu- 
minated section, the light appears plane polarized. If the 
model is now loaded and an external analyzer employed, the 
birefringence of the material between the illuminated sec- 
tion and the analyzer will produce interference effects 
that vary when either the model or the incident beam is 
moved about. If a section is illuminated at the surface 
near the observer and moved away from him, the space rate 
of formation of fringes, dN/ds, is* a measure of the bi- 
refringence in successive planes normal to the observation 
direction. Here H is the number of cycles of interfer- 
ence observed along the path and s is the path variable. 

Owing to the complications involved in accurately meas 
uring the movement of the light beam, another method of 
surveying the model is recommended. Let a plane polarized 
beam be projected into the model through a slit as previ- 
ously described. The scattering process in the model then 
acts as the analyzer rather than as the polarizer and no 
external analyzer is required. The value of dN/ds is 
now a measure of the birefringence in successive planes 
normal to the incident light beam. 

The relations between the birefringence and the stress 
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at any point in a loaded r.odel may "be considered as follows, 
When a material (such as Bakelite) "becomes optically aniso- 
tropic, its optical properties at a given point may "be con- 
veniently described in terms of three principal indices of 
refraction, n ly n 2 , and n 3 . These indices always lie at 
right angles to each other in space. If the optical ani- 
sotropy is due to stress, the three principal indices will 
lie in the direction of the three principal stresses (Si , 
So, and S 3 } ♦ Furthermore, the departures of the three 
indices n 2> and n 3 from the unstressed value n Q 

are proportional to 



S 1 , S 2 , and 


S 3 ; that i s , 


K(n 0 - n 1 ) = 


*» 


X(n 0 - n 2 ) = 




K(n 0 - n 3 ) = 




nt character 


istic of the 



In another form 

K(n 2 - n x ) = Si - So 

K(n 3 - n 2 ) = So - S 3 

K(n x - n 3 ) = S 3 - S x 

If S 1 <So<S 3 , the maximum shear is ;iven by 

K ( \ 
T max = g n n i - n 3 ' 

The maximum shear at a point may then "be found "by measuring 
tii " n 3 a ^ the point, provided that K is known. 

If a "beam of polarized light passes through the point 
in question al$ng the direction of S 2 and is polarized 
at 45° to S i and §« , the relative velocity of its com- 
ponents will depend on the ratio of gu to $u « Relative 
retardation will take place "between the components at a 
rate e;iven "by 



dR a (n 1 - n 3 ) 



or, expressing the relative retardation in wave lengths. 
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I ( ^ - = (S 3 - S x) ^3 - S x) 

where X is the wave length in air of the light used, 11 
is the amount of the relative retardation in wave lengths 
(cycles of interference), and C is another constant calle 
the stress-optical coefficient of the material . Inasmuch 
as it is proposed to analyze the light continuously along 
its path "by scattering, a fringe will appear each time the 
relative retardation equals a whole number of wave lengths; 
that is, when 

H = NA (N = 1, 2, 3, • . • m) 

If the direction of S 2 is unknown, the beam of light may 
"be passed through the point in various directions until 
dN/ds Is a maximum; that is, until the fringes appear most 
closely spaced, jDhe light beam is then traveling along the 
direction of S 2 . If it is assumed that the stresses re- 
main sensibly constant over a small increment of light path 
the formula may be written as 




If the increment is such that AIT * 1, then As is the 
fringe spacing and hence measurement of As determines 
(S 3 - S 1 ), If As is referred to as d, then 

Ss): - % = I 

Certain phases of the fringe formation process will 
now be considered. The scattered light is considered to 
be due to the vibration of submi cro s copi c particles in the 
scattering medium. These particles vibrate in the direc- 
tion of the electric vector of the incident beam of light. 
In birefringent material, such a particle vibrates under 
the influence of each of the two components into which the 
light divides. Light being a transverse wave motion, only 
the component of vibration normal to the direction of ob- 
servation yields visible light. Hence, the vibrating par- 
ticle should be viewed from a direction such that both 
light components take part in forcing the vibration. Inter 
ference between components will then be visible. If the 
light beam travels along S 2 and the observer looks along 
s i or s 3 1 n ° interference effects will be visible but, 
if the observer looks along a direction 45° to S 1 and 
S 3 , fringes with maximum visibility will be seen. The po- 
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sitions of minimum fringe intensity give the directions of 
5^ and S 3 . 

As a simple example of this procedure, consider the 
case of a tension member of rectangular cross section. 
Let a "beam of light "be collimated, plane polarized, and 
projected into such a model through a. slit. Let the slit 
lie parallel to the stress direction so that the illuminat- 
ed section is a longitudinal one through the piece. Let 
the light he polarized in a direction at 45° to the direc- 
tion of the stress. If the model is observed from a direc- 
tion either parallel to or normal to the direction of polar 
ization and at 90° to the incident light direction, fringes 
will appear which arc equally spaced. Hence d.Il/ds is a 
constant, and therefore the "birefringence is constant along 
the light paths. The result then indicates that the stress 
in the model is constant over the cross section. Such a 
fringe pattern is Shown in figure 1. A vector diagram of 
the light components is shown in figure 2. This diagram 
refers to the point at which light enters the model and "be- 
fore relative retardation has taken place between its com- 
ponents. Observation along A corresponds to the case in 
the plane-stress polariscope in which the polarizi?ig units 
are set for a dark field and, along 3, to the case of a 
light field. 

The principal stresses at a free "boundary may be indi- 
vidually measured since the stress perpendicular to such a 
"boundary vanishes and the other two stresses lie in the 
tangent plane. If the direction of the light is made par- 
allel to this tangent plane, the model may "be rotated about 
a point in the tangent plane until the maximum and the min- 
imum fringe spacings are obtained. These spacings will 
give the separate principal stress values. 



APPARATUS 



The apparatus employed for this study included a po- 
lariscope head, containing a light source, condensing 
lenses, a filter, a polarizing plate, and a slit. Above 
the head was placed a glass tank filled with a suitable 
liquid in which the models were immersed. A camera was 
provided for photographic recording. The arrangement of 
these elements is shown in figures 7* to 5. 



The light source for these studies must be intense 
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and a type H-6 Westinghouse 1 f OQO-watt mercury arc lamp of 
the water-jacketed variety was obtained. This lamp has 
been very satisfactory as regards intensity. Its spectrum, 
however, contains many more lines than the ordinary mercury 
arc and apparently has a fairly intense continuous back- 
ground , so that difficulty was experienced at first in ob- 
taining a monochromatic "beam. Such a "beam is highly de- 
sirable to prevent overlapping of various colors at high 
orders of interference, The combination of a TTratten No. 
77 filter with a Wratten 33 filter seemed best for photo- 
graphic recording and the Wratt en Bo. 77 alone for visual 
observation. The problem is further complicated by the 
variation in monochromatic scattering power for Bakelite, 
which is apparently strongest in the blue region, as would 
be expected. The lamp requires a separate transformer op- 
erating on 110 volts alternating current and a water sup- 
ply of about 1 gallon per minute. Provision for automatic 
shut-off of the power supply in the event of wat er- supply 
failure was not provided, although the automatic feature 
is very desirable. 

Two condensing lenses 6.5 inches in diameter served to 
collimate the light beam, the lamp being placed at the fo- 
cal point of the combi nat i on . The approximate speed of the 
condensing-iens system was f :G*83 « 

The polarizer employed was a sheet of type I Polaroid. 
No quarter-wave plate was used during the tests to be de- 
scribed, but one may occasionally be necessary to project 
circularly or elliptically polarized light into the model. 

The slit was approximately one-tenth inch wide and 
equal in length to the diameter of the collimated beam. 
Naturally, the intensity of the fringes and the speed of 
photographic recording vary directly with the slit width. 
The width used resulted in fringes easily visible in day- 
light. The exposures were of the order of 1 minute with 
the filter combination mentioned and a high-speed film. 
On the other hand, the slit width effectively determines 
the gage length over which the measurements are made; 
hence, it should be as narrow as possible consistent with 
satisfactory visibility and photographic exposure. Good 
visibility was obtained with the slit width greatly reduced, 
probably to about 0.01 inch. 

In order to avoid refraction at model surfaces, the 
model must be immersed in a clear liquid having the same 
refractive index (n = 1.57) as Bakelite. A mixture of 
the following proportions (by volume) 
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Halo wax oil No. 1000 (n = 1.6?) 



6 5 p e r c e n t 



Clear mineral oil (n = 1.43) 



34 percent 



was found to ^ive the correct refractive index. This liq- 
uid was contained in a gjass-'wall tank 3 inches by 3 inches 
by 5 inches deep, mad© of plate glass united with Sa.uer- 
ei sen cement . 



A beam in uniform "bending was selected as an example 
of the method. The stresses under a given loading were 
pho t o elast i cally measured and were subsequent ly checked by 
calculation according to the elementary theory. Two "beams 
were cut from a block of Eakelite BT-61-893. No annealin? 
was attempted. The models were 0.75 inch by 0.75 inch in 
cross section and were 6 inches in length. They were load- 
ed in a small fixture constructed for the purpose. Figures 
6 and 7 show a beam In place. Loads were found by measur- 
ing the deflection of the spring contained in the fixture. 

One beam was cut from the Bakelite block alon§ an. 
edge exposed during the "curing" process. This beam showed 
the prominent initial stress pattern seen in figure 3(a). 
In figures 8(b) and 8(c) are shown the fringe patterns in 
this beam due to one-half load and full load in bending. 
It will be noted that the distortion of the pattern from 
symmetry about the neutral axis decreases with increasing 
load, that is, the relative influence of the initial stress 
es decreases. 

Figure 9 shows the bending pattern of the second beam, 
which was cut from the parent material in such a way as to 
be well removed from previously exposed surfaces. This 
2^attern shows satisfactory symmetry and was taken as the 
basis for comparison with theory. Table I shows the values 
obtained. The procedure in obtaining the data was as fol- 
lows,: The negative of figure 0 was placed on the stage of 
a micrometer comparator and traverses were made along a 
vertical line at the center of the beam. Positions of the 
dark fringes are recorded in the first column of table I. 
Inasmuch as the image on the negative was reduced in size 
as compared with the beam itself, the readings were multi- 
plied by a factor (1.93) after subtraction from the initial 
reading as shown in the second column. The second column 
then gives the displacement of the fringes from the edge of 
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the "beam through which the light entered. The stress tralue 
obtained from the spacing "between a given pair of fringes 
was assumed to correspond to a point midway "between them* 
The third column gives the location of these points. The 
actual spacing, of the fringes is shown in the fourth col- 
umn. In the last column are recorded the stresses obtained 
by dividing the fringe value of 35 pounds per square inch 
per fringe by the spacing. 

The neutral axis of the beam (center line of the model 
represents a barrier across which fringe- spacing measure- 
ments cannot be made. Such a measurement would ^ive the 
spacing of a fringe with respect to itself and such a meas- 
urement has no meaning. The spacing across a region in 
which a reversal from tension to compression, or vice versa 
takes place bears no relation to the stress but dopends 
rather on the relative phase of the light components as 
they cross the neutral section. Only when this phase dif- 
ference is accidentally zero, will the measurement be sig- 
nificant in terms of stress. 

The values from table I are plotted in figure 10. 
Inasmuch as this feype of stress distribution is a straight 
line to a good approximation, such a line was drawn through 
the points. Calculation using the applied load of 350 
pounds and moment arms of 1.37 inches gives 

CT = § = * ^||— ; 6 1 * j£- = 3,400 pounds per square inch 

where cr is the maximum stress in the beam, M is the 
bending moment on the beam, and Z is the section modulus 
of the beam. This value compares with the value of 3 1 600 
pounds per square inch given by the intersection of the 
plotted value line with the model edges. The discrepancy 
here is 5.5 percent of the cornput ed value. The average 
deviation of the significant points from the curve is 4 
percent . 

As some of the equipment used was incapable of precise 
measurement, the result is considered good agreement. 
Sources of error may have been: 

(1) Error in model dimensions. 

(2) Error in spring calibration. 



(3) Error in spring-deflection measurement. 
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(4) Error in measuring roller positions on "beam. 

(5) Friction in loading fixture. 

(5) Error due to deflection of "beam. 

(7) Errors due to photographic distortion. 

(8) Errors in centering fringe line in micrometer 

comparat o r . 

(9) Error due to proximity of load points. 
(10) Error due to strains in tank "bottom* 
(11 ) Error due to strains in Polaroid disk. 

(12) Error due to lack of perfectly monochromatic 

light . 

(13) Initial optical ani sotropy in the model. 

(14) Error. in fringe value for Bakelite due to tem- 

perature . 

It is probable that items (3), (4), (5), (6), (9), 
and (13) contributed the larger part of the final error. 
An exception is the uncertainty in the reading of the po- 
sition of the first fringe with respect to the model edge 
This value was rejected as "being accounted for by items 
(10), (11), and (13), out chiefly "by (10). It: Is advi sa- 
ble to approach a 'bounding sur f ace from within and to re- 
duce the effect of these items "by plotting a curve such 
as has "been done here. Items (10) and (ll ) can affect on 
ly the position of the first fringe. 

In general, it may he said that the errors listed in 
items (l) through (9) are not characteristic of this meth 
od of stress analysis and may he controlled "by well-known 
methods. Errors due to (10) , (ll), and (12 ) can he elim- 
inated without difficulty "by the use of wellr.de signed 
equipment. Item (13) is controllable by suitable choico 
of model material; by careful machining; and by testing 
the model soon after cutting; or, when such testing is in 
convenient, by measuring the initial conditions without 
load and subtracting the result from later readings. The 
variation due to item (14; is small because of the manu- 
facturer ! s care in processing (reference 6). 
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An additional possi/ble source of error is due tc the 
assumption that the stress computed "by the spacing of two 
adjacent fringes exists at a point midway "between their:. 
This assumption is more nearly true as the spacing of 
fringes oecoir.es more uniform with high loads. In the pres 
ent example, the departure of the representative point 
from the midway position could have "been computed, "but co;. 
putation is generally impossible. Here, as in plane pho- 
toelastic analysis, a high fringe order is conducive to 
precise results. 

As an example of the usefulness of the method in a 
slightly more complex problem, the variation in the stress 
pattern due to a vertical hole through the center of a 
"beam is presented. A hole of 5/lS-inch diameter was 
drilled through the center of the beam previously analyzed 
Figures 11 and 12 show the stress patterns in the central 
plane and in the face of the beam, respectively. 

The disturbance in the fringe pattern at the load 
points was insignificant for both beams tested. This re- 
sult is in contrast to the prominent effects obtained in 
the usual j^lane- s t r e s s pattern. In the three-dimensional 
method, the components of stress that lie in horizontal 
planes are effective and no vertical stress components in- 
fluence the fringe positions. If the model had been il- 
luminated by a horizontal beam, the resulting fringes 
would have indicated the effect of the concentrations very 
prominently. The operator must not expect to see the 
types of fringe patterns that have come to bo associated 
with photoelastic methods in the past because the present 
procedure involves a new type of st r e's s-oi^t i cal behavior. 



CONCLUDE REMARKS 



The foregoing results indicate that the method devel- 
oped for three-dimensional photoelastic stress analysis 
is sound, that the apparatus described is suitable for a 
wide variety of stress investigations, and that the exper- 
imental techniques employed yielded satisfactory results. 
This method appears to be a substantial improvement over 
previous attempts to solve spatial stress problems by the 
use of photoelastic models and, in a large number of cases 
appears to be the only method thus far developed that will 
provide a solution in a convenient manner. 
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In the solution of the case of the "beam in simple 
"bending, as analyzed "by the method outlined in the paper, 
the agreement with accepted "beam theory was within 5,5 per-* 
cent, Ho basic assumptions are made in developing the the- 
ory other than those already well established in plane pho- 
toelastic procedure. 

Persons should familiarize themselves with the elemen- 
tary phenomena associated with the scattering of light "be- 
fore attempting experiments "by this method. 

It may at times he desirable to submit to >-ome of the 
disadvantages of the fixation method and, having "frozen 11 
a stress system into a model, to examine it by the methods 
described in this paper. The tension model used in this 
investigation was so treated. 

The principal stresses at a free boundary may be in- 
dividually measured since the stress perpendicular to such 
a boundary vanishes and the other two stresses lie in the 
tangent plane. If the direction of the light is made par- 
allel to this tangent plane, the model may be rotated about 
a point in the tangent plane until the maximum and the min- 
imum fringe spacings are obtained. These spacings will 
give the separate principal stress values. 

Many problems suggest themselves as suitable subjects 
for investigation. Certainly stress concentrations in 
three-dimensional systems should be thoroughly explored, 
particularly in connection with torsional stresses. It 
also seems possible that some of the problems which have 
been treated as two-dimensional ones in the past may not 
in fact be two-dimensional and that more accurate results 
could be obtained by the suggested method. 

Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va . , October 26, 1939, 
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TASLE I 



Data on Cross Section of Ba]:elit< 
BT-61-89.7 Beam, Correctly Cut 



Po si t i on of 
dark fringes 
on negative 
(in.) 



0.524 
.534 
.548 
.562 
.578 
.597 
.619 
.544 
.691 



edge 



. 758 
.796 
. 822 
. 843 
. 860 
.8 74 
. 890 
. 904 
.912 



Position 
of fringes 
in no del 
(in. ) 



0 

.019 
.046 
.073 
.104 
.141 
.183 
.232- 
. 322 
Gent e 
.45 2 
.5 25 
.575 
.616 
.548 
.676 
. 704 
. 733 
. 749 



Stress 
o cation 

(in. ) 


_ 

Fringe 
spa ci ng 

(in. ) 


Stress 
value s 

( lb . / sq . 


0. 010 


{ ) . uiy 


% , u JU 


. 032 


.02 7 


3 ,150 


.060 


.027 


3 ,150 


.089 


.031 


2 , 750 


.122 


.037 


2 , 300 


.15 2 


.042 


2,000 


.205 


.049 


1 ,750 


. 277 


.090 


950 



line of model 

.488 
.550 
• 5 9 5 
.5 32 
.552 
.690 
. 713 



.073 
.05 0 
.041 
.032 
.023 
.028 
.029 



1 , 150 

1 , 700 

2 ,050 

2 , 550 

3 ,05 0 
3,050 
2 ,950 




Figure 3.- General view of apparatus showing location of polariscope, tank 9 model, and camera. 
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Figure 4.- Internal layout of pol&ri scope. 
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Fig. 7 
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Fig. 8 

































(a) Initial stress pattern(no load).. 
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I Half bending load. 





































(c) Full bending load. 



Figure 8.* Stress patterns in a Bakelite beam showing initial stress 
resulting from an edge that has been exposed during curing. 
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Figure 11.- Beam with vertical hole, center section stress pattern. 




Figure 12.- Beam with vertical hole, edge section stress pattern. 



